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Abstract. We consider a matrix pencil whose coefficients depend on a posi- 
tive parameter e, and have asymptotic equivalents of the form ae A when e goes 
to zero, where the leading coefficient a is complex, and the leading exponent 
A is real. We show that the asymptotic equivalent of every eigenvalue of the 
pencil can be determined generically from the asymptotic equivalents of the 
coefficients of the pencil. The generic leading exponents of the eigenvalues are 
the "eigenvalues" of a min-plus matrix pencil. The leading coefficients of the 
eigenvalues are the eigenvalues of auxiliary matrix pencils, constructed from 
certain optimal assignment problems. 



Perturbation de valeurs propres de faisceaux matriciels 
et probleme d'affectation optimale 

Resume. Nous considerons un faisceau matriciel dont les coefficients, dependant 
d'un parametre e, ont des equivalents asymptotiqucs de la forme ae A , lorsque 
e tend vers zero par valeurs positives, le coefficient dominant a etant complexe 
et l'exposant dominant A etant reel. Nous montrons qu'un equivalent asymp- 
totique pour chacune des valeurs propres du faisceau peut etre determine 
generiquement a partir des equivalents des coefficients du faisceau. Les ex- 
posants dominants des valeurs propres sont les valeurs propres d'un faisceau 
matriciel min-plus, et les coefficients dominants sont les valeurs propres de fais- 
ceaux auxiliaires, construits au moyen de problemes d'affectation optimale. 



Abridged French version 

Nous considerons un faisceau matriciel A e = A e ,o + XA e ,i + • ■ ■ + X d A t) d, oil 
X est une indeterminee, et oil pour tout < k < d, A e>k est une matrice n x 
n dont les coefficients, {A e> h)iji sont des fonctions continues a valeurs complexes 
d'un parametre positif e. On s'interesse au comportement asymptotique, lorsque e 
tend vers 0, des valeurs propres C e de A e , qui sont par definition les racines du 
polynome det(.A £ ). Nous supposons que pour tout < k < d, on a des matrices 
a k - ((a k )ij) G C nx " et A k = {{A k ) tJ ) e (K U {+oo}) nx ™ telles que {A e<k )ij = 
(ak)ij^ Ak ^ ij + o(e^ Ak ^ ij ), pour 1 < i,j < n. Lorsque (A k )ij = +oo, ccla signifie, par 
convention, que {A t ,k)ij est nulle au voisinage de 0. Nous cherchons un equivalent 
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asymptotique de la forme C e ~ Ae A , avec A G C \ {0} et A £ R, pour chaque valeur 
propre C e de A e . 

Lorsque A e = A e ,o ~ Xid, ou id designe la matrice identite, il s'agit du probleme 
classique de perturbation de valeur propres de matrices. En particulier, lorsque la 
matrice A e ,o est de la forme A e ,o = Ao t o + eb, pour une matrice complexe b = (bij), 
une theorie initiee par Visik et Ljusternik VL60 , et completee par Lidskii Lid65^ 
identifie les exposants dominants A, pour des valeurs generiques des parametres 
bij, et donne des formules explicites pour les coefficients dominants A, faisant inter- 
venir les valeurs propres de certains complements de Schur construits a partir de 
b. Voir [MBQ97] pour une presentation recente. Cependant, pour certaines valeurs 
de b, ces complements de Schur peuvent ne pas exister. De tels cas singuliers ont 
ete etudies en particulier dans |MBQ97I lME98| . Dans |ABG04j . nous avons domic 
une premiere extension du theoreme de Lidskii |Lid65| , donnant sous certaines con- 
ditions structurelles les equivalents des valeurs propres de la matrice A e ,o- Nous 
etendons ici ce resultat au cas des faisceaux, ce qui permet de resoudre les cas qui 
demeuraient singuliers dans ABG04 . 

Afin d'enoncer le resultat principal, rappelons quelques notions d'algebre min- 
plus. Le semi-anneau min-plus, R m i n , est l'ensemble R U {+00} muni de l'addition 
(a, b) <— > a b :— min(a, b) et de la multiplication (a, b) 1— > a (g> b := a + b. On ecrira 
parfois ab au lieu de a eg) b. On notera (D := +00 et 1 := le zero et l'unite de R m i n) 
respectivement . 

Nous associons au faisceau A e le faisceau matriciel A = Aq © XA\ © ■ • ■ © X d Ad 
a coefficients dans M m i n - Les coefficients de A, Aij, sont done des polynomes 
formels a coefficients dans K m i n , en l'indeterminee X. Nous appelons polynome 
caracteristique min-plus le permanent Pa = perm A. Rappelons que pour une ma- 
trice B = (B^) de taille n x n, a coefficients dans un semi-anneau quelconque, 
le permanent de B est defini comme la somme sur toutcs les permutations a du 
poids \o\b = ^io-(i) ' ' ' B n(J („). Lorsque les coefficients de B appartiennent a M. m i n , 
\o\b = B la (i} + • • • + B na ( n j, et permi? est la valeur d'une affectation optimalc 
dans le graphe value associe a B. 

Si P est un polynome formel a coefficients dans M m i n , on note P la fonction 
polynome associee a P. L'application P 1— > P est un cas particulier de transformee 
de Legendre-Fenchel |BCOQ92| § 3.3.1]. Cuninghame-Green et Meijer |CGM80| ont 
montre que la fonction polynome P(x) peut se factoriser de maniere unique sous la 
forme P(x) = a(xffici) • ■ • (x®c n ), avec a, ci, . . . , c n £ K m ; n . Les nombres c\,... 1 c n 
sont appeles racines de P. Nous noterons 71,... ,7jv les racines du polynome car- 
acteristique Pa, aussi appelees valeurs propres du faisceau A. Les racines ji, et done 
la fonction polynome Pa, peuvent etre calcules en temps 0(n A d) en s'inspirant de 
la methode de Burkard et Butkovic BB03 . 

Pour toute matrice B £ telle que permi? ^ (D, nous definissons le graphe 

Opt(i3) forme des arcs participant a une affectation optimale : les noeuds de Opt(i?) 
sont 1, . . . ,n, et il y a un arc de i a j s'il y a une permutation a telle que j = a(i) 
et \<t\b = perm B. 

Nous dirons que deux vecteurs U, V de dimension n a coefficients dans R m i n \ {©} 
forment une paire hongroise relativement a B si, quels que soient i,j, on a B^ > 
UiVj, et si U\ ■ ■ ■ U n Vi ■ ■ ■ V n = permB, les produits s'entendant dans le semi- 
anneau min-plus. Ainsi, (U, V) n'est autre qu'une solution optimale du probleme 
lineaire dual du probleme d'affectation (see |BR97I 2.4]). En particulier, une paire 
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hongroise existe des que perm B ^ (D, et elle peut etre trouvee en temps 0(n 3 ) grace 
a l'algorithme hongrois. Pour toute paire hongroise (U, V), on definit le graphe de 
saturation, Sat(£>, U, V), qui a pour noeuds 1, . . . ,n, et pour arcs les couples 
tels que B^ = UiVj. 

Pour chaque racine finie 7 de Pa, nous definissons les graphes 
Opt (7), . . . , Opt d (7) : Opt fe (7) a pour nceuds 1, ...,n, et a un arc de i a 
j si £ Opt(A(7)) et (Ak)ijj k = Ay (7). Pour toute paire hongroise 

(U, V) relativement a la matrice ^(7), nous definissons aussi les graphes 
Sato(7, U, V), . . . , Satd(7, U, V), qui sont obtenus en remplacant Opt(A(^)) par 
Sat(A(7), U, V) dans la definition de Opt (7), . . . , Opt d (7). Enfin, si G est un 
graphe ayant pour noeuds 1, . . . , n, et si b 6 C nxn , nous definissons la matrice b G , 
telle que (b G )ij — bij si s G, et (6 )ij = sinon. Dans 1'enonce du theoreme 
qui suit, les valeurs propres sont comptees avec leurs multiplicites. 

Theoreme 1. Soit 7 une racine finie du polyndme caracteristique min-plus Pa- 
Pour chaque < k < d, notons G k le graphe egal a Opt fe (7) ou bien a Sat/-(7, U, V), 
pour un choix quelconque de la paire hongroise U,V relativement a ^(7). Con- 
siderons le faisceau auxiliaire := a °+Xa 1 1 + - ■ -+X d a d d . Alors, si le faisceau 
a m 7 valeurs propres non-nulles, Ai, . . . , A m ^ , le faisceau A € admet m 7 valeurs 
propres jC £i i, . . . ,C e ^ m ^ ayant des equivalents respectifs de la forme C tt i ~ A^e 7 . En 
outre, si est une valeur propre de multiplicity m' du faisceau le faisceau A e 
admet m 7 valeurs propres supplementaires C e telles que e _7 £ e converge vers 0, et 
toutes les autres valeurs propres C e de A e sont telles que le module de e~ 7 £ e tend 
vers Vinfini. Enfin, pour des valeurs generiques des parametres (afe)y, to 7 coincide 
avec la multiplicity de la racine 7, et m 7 coincide avec la somme des multiplicites 
des racines de Pa strictement superieures a 7. 

Dans la version en anglais de la presente note, nous illustrons le Theoreme \T\cn 
raffinant des resultats de Najman [Naj99 . 



We consider a matrix pencil of the form 

A £ = A e<0 + XA e s + ■■■ + X d A t4 , 

where, for every < k < d, A e< k is a n x n matrix whose coefficients, (A et k)ij, are 
complex valued continuous functions of a nonnegative parameter e, and X is an 
indeterminate. We are interested in the asymptotic behavior, when e tends to 0, of 
the eigenvalues C e of A e , which are the roots of the polynomial det(^l e ). 

We shall assume that for every < k < d, matrices ak = ((afe)y ) G C nx ™ and 
A k = {(Ak)ij) 6 (KU {+oo}) nx ™ are given, so that 

{A,k)a = (akhe {Ak) ^ + o(e^ fe )«) , for all 1 < i,j < n , 

when e tends to 0. When (Ak)ij — +00, this means by convention that {A e< k)ij is 
zero in a neighborhood of zero. We look for an asymptotic equivalent of the form 
£ £ ~ Ae A , with A e C \ {0} and AgK, for every eigenvalue £ e of A e . 

When A e = A e .Q — Xid, where id denotes the identity matrix, we recover the 
classical problem of perturbation of eigenvalues of matrices. In particular, when 
A e ,o is affine in e, i.e. A e ,o = Aq,q + eb, for some complex matrix b — {bij), a theory 
initiated by Visik and Ljusternik |VL60| and completed by Lidskii |Lid65j identifies 
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the leading exponents A for generic values of the parameters bij , and gives explicit 
formulae for the leading coefficients A, involving the eigenvalues of certain Schur 
complements constructed from b. See MB097 for a recent presentation. For some 
values of b, however, these Schur complements may not be defined. Such singular 
cases have received much attention, see in particular |MB097l ^IE98 . In ABG04 , 
we generalized the theorem of |Lid65| to the case where A t fl is not affine in e. We 
showed that the leading exponents A of the eigenvalues can be computed, gener- 
ically, as the "roots" of a certain min-plus polynomial constructed from A e . We 
also showed that under structural conditions (when certain bipartite graphs con- 
structed from the leading exponents of the entries of A €j o have a perfect matching), 
the leading exponents of all the eigenvalues can be obtained by Schur complement 
formulas which extend the ones of |Lid65| . However, when these structural condi- 
tions do not hold, the question remained open. We show in Theorem ^ below that 
the leading coefficients of the eigenvalues of the pencil A e are obtained, generically, 
as the eigenvalues of certain auxiliary pencils (independent of e), constructed by 
solving certain optimal assignment problems. This gives a complete generalization 
of the theorem of |Lid65| . We shall illustrate Theorem^n Corollary^ by refining 
results of Najman |Naj99| concerning a singular perturbation of an affine pencil. 

In order to state the main result, we need some min-plus algebraic constructions. 
Recall that the min-plus semiring, R m i n , is the set KU {+00} equipped with the 
addition (a, b) 1— > min(a, b) and the multiplication (a, b) 1— > a + b. We denote by "©" 
the min-plus addition, and by "<g>" or concatenation the min-plus multiplication. 
We denote by = +00 and 1 = the zero and unit elements of K m in, respectively. 
The max-plus semiring, R max , is obtained by replacing "min" by "max" and +00 
by —00 in this definition. 

We associate to the matrix pencil A e the matrix pencil with coefficients in M m m, 
A = Aq © XA\ © • ■ • © X d Ad- Here, the entries of A, Aij, are formal polynomials 
in the indeterminate X with coefficients in M m i n - We consider the formal min-plus 
characteristic polynomial Pa = perm A. Recall that for a n x n matrix B = (Bij) 
with entries in any semiring, the permanent of B is defined as the sum over all 
permutations a of the weight \<j\b = ^io-(i) ' ' ' Bna(n)- Thus, if B is a matrix with 
entries in R m i n , the weight of a, \<j\b, is the usual sum B 1<7 ^ + • • • + B n<T („), and 
perm B is the value of an optimal assignment in the weighted graph associated to 
B. (See |BR97| for more background.) 

If P is a formal polynomial with coefficients in R m i n , we denote by P the polyno- 
mial function associated to P. The map P 1— > P is a specialization of the Legendre- 
Fenchel transform |BGOQ92[ § 3.3.1]. Cuninghame- Green and Meijer |CGM80| 
have shown that the min-plus polynomial function P(x), can be factored uniquely 
as P(x) = a(x © ci) • • • (x © c n ). with a, ci, . . . , c n G R m in, where n is equal to 
the degree of P, degP. The numbers ci, . . . , c„, are called the corners of P. They 
coincide with the points of nondifferentiability of P. If c is a corner, the multiplicity 
of c, which is equal to the number of indices i for which c$ = c, coincides, when 
c 7^ D, with the variation of slope of P at c, P'(c~) — P'(c + ). The multiplicity of 
the corner (D is equal to the valuation of P, valP. We denote by 71,... ,7iv the 
corners of the characteristic polynomial Pa, that we call the (algebraic) eigenvalues 
of A. Note that the valuation val Pa can be computed by introducing the matrix 
valA G R„J n n , such that (val^4)y = val Ay. Then, valP^ is equal to the min-plus 
permanent of the matrix val A. By symmetry, the degree deg Pa is equal to the 
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max-plus permanent of the matrix deg^4 G such that (deg A)tj = degAjj. 

The corners 7j (and so, the polynomial function Pa) can be computed in (9(n 4 c?) 
time by adapting the method of Burkard and Butkovic BB03 . (It is not known 
whether the sequence of coefficients of the formal polynomial Pa can be computed 
in polynomial time.) 

For any matrix B 6 R^in™ sucn that permi? ^ (D, we define the graph Opt(B) as 
the set of arcs belonging to optimal assignments: the nodes of Opt(B) are 1, . . . ,n 
and there is an arc from i to j if there is a permutation a such that j = a(i) and 
\<j\b = permB. 

We shall say that two vectors U, V of dimension n with entries in R m i„ \ {©} 
form a Hungarian pair with respect to B if, for all i,j, we have Bij > UiVj, and 
Ui • • ■ U n Vi ■ ■ ■ V n = permi?, the products being understood in the min-plus sense. 
Thus, (U, V) coincides with the optimal dual variable in the linear programming 
formulation of the optimal assignment problem (see |BR97I 2.4]). In particular, a 
Hungarian pair always exists if the optimal assignment problem is feasible, i.e., if 
permi? ^ C, and it can be computed in 0(n 3 ) time by the Hungarian algorithm. 
For any Hungarian pair (U, V), we now define the saturation graph, S&t(B,U,V), 
which has nodes 1, . . . ,n and an arc from i to j if B^ = UiVj. 

For every finite corner 7 of Pa, we define the digraphs Opt (7), . . . ,Opt d (7): 
Opt fc (7) has nodes l,...,n, and an arc from i to j if (i,j) £ Opt(A(j)) and 
{Ak)ijl k — Aij(~f). For every Hungarian pair (U, V) with respect to the ma- 
trix ^(7), we also define the digraphs Sat (7, 17, V), ... , Satd(7, U, V) by replacing 
Opt(i(7)) by Sat(i(7), U, V), in the definition of Opt ( 7 ), . . . , Opt d ( 7 ). Finally, 
if G is any digraph with nodes 1,. . . ,n, and if b G C" x ", we define the matrix 
b G , which is such that (b )ij = bij if (i,j) 6 G, and (b G )ij = otherwise. In the 
following theorem, and in the sequel, eigenvalues are counted with multiplicities. 

Theorem 1. Let 7 denote any finite corner of the min-plus characteristic polyno- 
mial Pa- For every < k < d, let Gk be equal either to Opt fc (7) or Satfc(7, U, V), 
for any choice of the Hungarian pair U, V with respect to ^(7). Consider the aux- 
iliary pencil 

a (7) := a Ga + Xaf 1 +■■■ + X d a Gd . 

Then, if the pencil has m 7 non-zero eigenvalues, Ai , . . . , , the pencil A e has 
m 1 eigenvalues C e ,i, . . . , with respective equivalents of the form C e ^ ~ A^e 7 ; 

ifO is an eigenvalue of multiplicity m'^ of the pencil oW, the pencil A e has precisely 
m'^ eigenvalues C e such that e~ 7 £ £ converges to zero, and all the other eigenvalues 
C e of A e are such that the modulus of e _7 £ £ converges to infinity. Moreover, for 
generic values of the parameters (ak)ij, m 7 coincides with the multiplicity of the 
corner 7, and m 7 coincides with the sum of multiplicities of all the corners greater 
than 7. 

Since the sum of the multiplicities of the corners of Pa is equal generically to 
the degree of det(.4 c ), and since A e has a number of identically zero eigenvalues 
generically equal to the multiplicity of as a corner of Pa , Theorem ^ provides 
generically asymptotic equivalents of all the eigenvalues of A t . 
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Example 1. Consider A e = A e ,o — Xid, and 
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The associated min-plus matrix pencil and characteristic polynomial function are 
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where hj G C. 
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P A (x) = (i®0) 2 (iffll) 



so that the corners are 71 =72 = 0, with multiplicity 2, and 73 = 1, with mul- 
tiplicity 1. We first consider the corner 7 = 0. Then U — V — (0,0,0) yields a 
Hungarian pair with respect to the matrix 

" O01 Oo 
i(0)= Oo 0i 1 

. Oo 10! 

where we adopt the following convention to visualize the digraphs Satfc(i7, V): an 
arc belongs to Satfc(£7, V) if k is put as a subscript of the entry Aij(Q). For 
instance, Au(0) = 0, and (1, 1) belongs both to Sato(f7, V) and Sati(f7, V). Entries 
without subscripts, like ^23(0) = 1, correspond to arcs which do not belong to 
Sat([/, V). The eigenvalues of the auxiliary pencil cS ' are the roots of 
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= A(— A + Xbn + foi2&2i 

Theorem ^ predicts that this equation has, for generic values of the parameters bij, 
two non-zero roots, Ai, A2, which yields two eigenvalues of A e , C e ,m ~ A m e° = A m , 
for m = 1,2. Consider finally the corner 7 = 1. We can take U = (0,1,1), 
V = (—1,0,0), and the previous computations become 

i(l) = 

The latest equation yields A(6i 2 &2i + &13&31) + &12&23&31 + &13&32&21 - &21&12&33 - 
&31&13&22 = 0. Theorem n predicts that this equation has generically a unique 
nonzero root, Ai, and that there is a branch ~ Aie. 

As a typical application of Theorem ^ let us consider the following singular 
perturbation of an affine pencil, A £ = eX 2 m+Xc+k : already considered in |Naj99 . 
For non-zero values of the entries of the matrices m, c, and k, the associated min- 
plus characteristic polynomial function is Pa{x) = (0 © x)™(0 © lx) n . Moreover, 
the pencils Xc + k and Xm + c generically both have n finite non-zero eigenvalues, 
denoted by Ai, . . . , A n , and /ii, . . . , /i„, respectively. Then, it is easy to derive from 
Theorem n that the pencil A e has n eigenvalues C €l i ~ Aje , and n eigenvalues 



Consider now the following non-generic situation. Assume that the 



pencil cX + k is given, that it is regular, and that its Weierstrass normal form 
comprises qo Jordan blocks for the eigenvalue 0, with respective sizes sj,. 



and q a 
set do 



Jordan blocks for the eigenvalue 00, with respective sizes s^. 



°o ' 
We 



'0 > 



We also denote by q' the number of 
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one dimensional Jordan blocks for the eigenvalue of the pencil cX + k. We denote 
by Ai, . . . , X r the finite non-zero eigenvalues of cX + k (of course, r + do + doc = n). 
We also denote by fii, . . . , fit the finite non-zero eigenvalues of the pencil Xm + c. 
We say that an eigenvalue C e is of order e A if £ e ~ Ae A , for some A 6 C \ {0}. 
The following result should be compared with |Naj99| , where partial results are 
obtained in a similar situation. 

Corollary 1. The pencil A e = eX 2 m + Xc + k has precisely 

(i) r eigenvalues of order e , which converge respectively to \, for i = 1, . . . , r; 

(ii) t eigenvalues or order e _1 , which are respectively equivalent to \i{e~ x , for 
i = l,...,t. 

It has at least 

(hi) 2qo — q' eigenvalues identically equal to zero. 
Finally, for generic values of the parameters m^, we have t = n — q^, and the 
pencil A e has precisely: 

(iv) + 1 eigenvalues of order e~ 1 A s ot+ 1 ) ; for i = 1, . . . , q^; 

(v) Sq — 2 eigenvalues of order e 1 /( s o _2 ) ; for every i such that 1 < i < r and 
s >2. 

Corollary provides, for generic values of to, the leading exponents of all the 
eigenvalues of the pencil A e . In cases Irvkvl the generic values of the leading coef- 
ficients of the eigenvalues can be determined by formulae essentially similar to the 
case of |Lid65 AB G04| . This will be detailed elsewhere. 

Acknowledgement. The third author thanks Jean-Jacques Loiseau for having 
suggested to look for a generalization of the result of |AB(t04) to matrix pencils. 
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